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Neutrino Flux Bounds and Prospects for High Energy and Ultrahigh 
Energy Neutrino Source Detection 
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Abstract. After briefly reviewing various hadronic neutrino 
source models, we show how to construct generic upper flux 
bounds. We then turn to the problem of neutrino propaga- 
tion through the inner Earth and neutrino detection in water- 
based Cerenkov detectors. Applying the formalism thus de- 
veloped to the Mannheim-Protheroe-Rachen (MPR) and the 
Waxman & Bahcall flux bounds, we find that event rates of 
several hundred to thousand events per year might be pos- 
sible in next-generation neutrino telescopes. However, a to- 
mography of the inner Earth will face severe constraints due 
to the statistical error of the event rates to be expected. 



1 Hadronic neutrino source models 

Models for sources producing ultrahigh energy neutrinos can 
be divided into hadronic and non-hadronic ones, depending 
on whether protons play a key part in the neutrino produc- 
tion. In this work we shall limit ourselves to the former class 
of models, referring the reader to Learned and Mannheim 
(2000) for a discussion of non-hadronic source models such 
as the decay of relic particles. 

Independent of the specific source under consideration, 
the main idea behind hadronic neutrino production is that a 
beam of accelerated protons (taken to be of power-law form, 
rip(Ep) oc E~ s ) hitting a target produces pions, which sub- 
sequently decay and thus yield neutrinos. Here, the target 
may consist of either protons or photons. 

In case of a target made of protons, if the pions decay prior 
to any interaction with the surrounding medium, the resulting 
neutrino production spectrum is of the same form as the orig- 
inal spectrum of accelerated protons, i.e. QfF(E u ) oc E~ s . 
Otherwise, the neutrino spectrum is steepened (Rachen and 
Meszaros, 1998). 

An important example of a photon target is given by the 
cosmic background radiation, which leads to the Greisen- 
Zatsepin-Kuzmm cutoff. In neutrino sources, target photons 
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may be produced by means of synchrotron radiation. As- 
suming that the synchrotron photon spectrum is of the form 
N^(E 7 ) oc E~ a , one may show that the corresponding neu- 
trino spectrum is also a power law, QP~<(E V ) oc E v ^ °^ 
(Learned and Mannheim, 2000). 

In the following, we mention several examples of source 
models which have been proposed for ultrahigh energy neu- 
trino production: 

Terrestrial atmosphere: 

Neutrinos produced by cosmic -ray interactions in the ter- 
restrial atmosphere can be detected with present-day water- 
based Cerenkov telescopes and thus serve as an important 
benchmark (Domogatsky, 2001; Andres et al., 2001). 

Sun: 

The same processes responsible for neutrino production in 
the terrestrial atmosphere can also take place in the Sun. 
However, the neutrino flux thus to be expected is so low that 
it will only be discernible with next-generation neutrino de- 
tectors (Moskalenko et al., 1991; Hettlage et al., 2000). 

Crab nebula: 

Neutrons resulting from the photodisintegration of iron nu- 
clei may be accelerated in the outer gap of the Crab pulsar 
and may subsequently decay into protons. If these interact 
with matter in the nebula, pions and thus neutrinos can be 
produced (Bednarek and Protheroe, 1997). 

Supernova remnants: 

Protons may be shock accelerated and by means of colli- 
sions in the remnant may produce pions and thus neutrinos 
(Gaisser et al., 1998). 
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AGNs: 



In active galactic nuclei (AGNs), protons may be shock ac- 
celerated either in the accretion disk of the central black hole 
or in the jets. Depending on the site of acceleration, either 
the thermal photons of the accretion disk or the synchrotron 
radiation in the jet may serve as the target needed for pion 
production (e.g. Nellen et al. (1993); Stecker and Salamon 
(1996)). 

It should be noted that the detection of neutrinos emitted 
by AGNs would be a smoking gun for hadron acceleration. 

Galaxy clusters: 

The ionized gas in a galaxy cluster might act as the target for 
cosmic rays escaping from cluster galaxies (Colafrancesco 
and Blasi, 1998). 

GRBs: 

Protons might be accelerated in gamma-ray bursts (GRBs) 
(e.g. Waxman and Bahcall (1997, 2000); De Paolis et al. 
(2001)). The main advantage of gamma-ray bursts is that 
they allow a precise timing, which reduces the background 
considerably. In addition, assuming that two neutrinos emit- 
ted at the same moment in a GRB located at a distance d 
arrive at different times t and t + At, one may compute the 
speed difference Av u of the two neutrinos by means of 
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Hence GRBs could offer a possibility for checking neutrino 
dispersion. 

Fig. 1 shows the muon neutrino fluxes predicted for vari- 
ous source models, if no neutrino flavor oscillations are as- 
sumed. 



2 Neutrino flux bounds 

The multitude (and uncertainty) of neutrino source models 
suggests the question whether one can obtain a generic up- 
per neutrino flux bound, thus constraining the number of 
events to be expected from future neutrino telescopes. In 
order to construct such a bound, we start by assuming that in 
a neutrino source all protons are confined within the accel- 
eration region (while neutrons and neutrinos may escape), 
that there is a power law target photon spectrum with in- 
dex a, n 7 (.E 7 ) oc E~ a , and that the upper limit on the 
cosmic ray proton spectrum as observed on Earth has the 
form n Pt0hs (E p ) oc E~ 2 - 75 between 3 x 10 6 and 10 12 GeV 
(Mannheim et al., 2001). 

Particle physics considerations then show that the neutrino 
production rate Q v within the source is given by 
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Fig. 1. Muon neutrino and antineutrino fluxes for 3C 273 (Nellen 
et al., 1993; Learned and Mannheim, 2000), the Crab nebula 
(Bednarek and Protheroe, 1997), the supernova remnant SNR IC 
443 (Gaisser et al., 1998; Learned and Mannheim, 2000), and the 
Coma cluster (Colafrancesco and Blasi, 1998). No neutrino flavor 
oscillations are assumed. 
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Fig. 2. Generic muon neutrino and antineutrino flux bound for 
sources with low ( r n7 < 1) and high neutron opacity ( r„ 7 ^> 1 ) 
(MPR flux bounds), and generic bound for neutrinos emitted by 
quasars (quasar), as obtained by Mannheim et al. (2001). The 
shaded region is the allowed range for the MPR flux bounds. Also 
shown are the generic GRB flux bound given by Mannheim (2001) 
(GRB) and the flux bound obtained by Waxman and Bahcall (1999) 
(WB). No neutrino flavor oscillations are assumed. 
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where Q n denotes the neutron production rate (Mannheim et 
al., 2001). The source contributes to the cosmic ray proton 
flux, if the neutron escape probability P eS c,n does not vanish. 
Indeed, the cosmic ray production rate can be estimated as 
(Mannheim et al., 2001) 



Qcr{Ep) ~ Qn(Ep) X P csc ,n(Ep) 



Qn(E p ) X 



. 1 _ e " 
T 

(T< 1) 
(T>1) ' 



-Qn(Ep) 



(2) 



Here r constitutes the neutron optical depth for the source 
radius. Finally, we note that the bolometric photon and neu- 
trino luminosities are related by L 7 = 2L U for a = 1 and 
L 1 « L v for a = (Rachen and Meszaros, 1998). 

We can now give a formal recipe for constructing generic 
upper flux bounds (Mannheim et al., 2001): 

1 . Get the neutron production rate Q n . Note that due to the 
rescaling in step 4, Q n need to be known up to a factor 
of proportionality only. 

2. By means of Eqs. 1 and 2 compute the corresponding 
neutrino and cosmic ray proton production rates (i.e. 
Q v » and Q cr ). 

3. Integrate the neutrino and cosmic ray proton production 
rates over all sources in the universe to obtain the neu- 
trino and cosmic ray proton intensities I at the Earth. 
Formally, this may be done using the formula 



M(E, z) 



(1 + zf dV c dP s( 



4< 



dz dV c 
Q((l + z)E,z)dz, 



where z denotes the redshift, the luminosity dis- 
tance, V c the comoving volume, and dP sourco /dl4 the 
source density, which analogously to Q n in step 1 needs 
to be known only up to a factor of proportionality. The 
modification factor M takes into account the attenua- 
tion processes along the line of flight. For neutrinos, 
M(E 7 z) — 1. The factor 1 + z in the argument of Q 
describes the energy degradation due to the expansion 
of the universe. 

4. Normalize 7 cr and thus I v so that 7 cr is tangent to the 
observational upper limit of the cosmic ray proton inten- 
sity or so that 7-rays aren't overproduced (whichever is 
the more stringent condition). 

It should be noted that the bounds thus computed pertain to 
hadronic source models only, so that there might be an ad- 
ditional contribution from non-hadronic sources. One may 
now use this recipe for computing various bounds: 

Generic bounds (Mannheim et al., 2001) are obtained, if 
one assumes that the neutron production spectrum is given 
by Q n (E n ) oc E^ 1 cxp(-E n /E max ) with a cutoff energy 
B max and that the source distribution follows that of galax- 
ies and AGNs. The corresponding cosmic ray proton flux 



as a function of E max is fitted to the observational bound. 
We shall refer to these bounds as the Mannheim-Protheroe- 
Rachen (MPR) flux bounds. 

In order to get a bound for EGRET blazars (Mannheim et 
al., 2001), we assume that the production rate of accelerated 
protons is of the form Q p (E p ) oc E~ 2 cxp(— E p /E max ), 
corresponding to non-relativistic Fermi acceleration, and that 
the spectral index of the target photon spectrum has the value 
a = 1. From Q p and the proton loss time scale one may com- 
pute the proton spectrum N p in the source, and from N p and 
the timescale of pj — ► n the neutron production rate Q n , 
which may serve as the starting point for the general recipe 
outlined above. Here, as the production rates depend on the 
luminosity, the average over the luminosity function dN/dL 
is used: 



(Q^,cr(E)) 



jQ^ CI { El L)^dL 



For GRBs (Mannheim, 2001) one takes the spectrum of 
accelerated protons to be of the form Q P (E P ) oc E p 2 , but 
assumes that the photon target spectral index is given by a — 
0. The calculation can then be carried out analogously to that 
for the blazar bound. 

Finally, if it is assumed that the proton input spectrum is of 
the form Q p (E p ) oc E~ 2 and that the sources are optically 
thin for neutrons, one obtains the Waxman & Bahcall bound 
(Waxman and Bahcall, 1999). 

In Fig. 2 various upper flux bounds for the muon neutrino 
flux from hadronic sources are shown. 



3 Event rates in neutrino telescopes 

The decay of a charged pion leads to the production of one 
electron and two muon neutrinos, where here and in the fol- 
lowing no distinction between particles and antiparticles is 
made. This implies that initially the ratio of the various neu- 
trinos should be given by Q Ve : Q Vfx : Q Vr =1:2:0. 

However, during the flight to Earth neutrino flavor oscil- 
lations should occur, leading to an equipartition among the 
various flavors (Learned and Pakvasa, 1995; Husain, 2000). 
Hence the flux ratio as measured on Earth should have the 
value 

I„. : 7^ : I Vt = 1 : 1 : 1, 

so that a considerable amount of tauon neutrinos is to be ex- 
pected, although virtually no tauon neutrinos are produced 
within neutrino sources. Neutrinos reaching a detector from 
below have to cross the inner Earth first. This propagation 
can formally be described by means of the set (Hettlage and 
Mannheim, 2001) 



dl vi 
dt 



-Cj,tot 



da k ^(E',E) 
dE 



Iu k (E% 0) 



of cascade equations, where ui denotes the total neutrino 
cross section of flavor i, the sum runs over all flavors, and 



4 



T n7 >>l 

T„<1 



0.4 0.6 
cos(nadir angle) bins 




0.4 0.6 
cos(nadir angle) bins 



Fig. 3. Event rates of muon (solid) and tauon neutrinos (clashed lines) for bins of the cosine of the nadir angle. Left: Event rates for the MPR 
flux bounds with high ( r„ 7 ^> 1) and low neutron source opacity ( r„ 7 < 1). Right: Event rates for the Waxman & Bahcall flux bound. 



ddk^i /dE covers both neutrino-nucleon and leptonic neu- 
trino interactions, and the decay of tauons. For the parton 
distributions we use the CTEQ5DIS functions (Lai, 2000), 
and we discretize the transport equations. The resulting set 
of equations can then be solved semi-analytically. 

As a rough approximation to the event rate N due to neu- 
trino detection in water-based Cerenkov telescopes, we em- 
ploy the formula 



Ni(E) 



pA, 



cir 



d(T CC 

dE'^ r I i {E')L det 



E 



?l da 



cc 



dE' 



+ J^xJ^ dE'^(E^(E',x))I t (E^(E',x))j.(4) 

Here again i denotes the flavor. A cS constitutes the effective 
detector area, L the detector height, p the matter density, and 
m p the proton mass. E^ (E 1 , x) is the energy (of the lepton 
i) which is degraded to E' along a path of length x. In the 
following, we shall assume A c k = 1 km 2 , L = 1 km, and 
p = 1 g/cm 3 . 

Whereas the first term on the right hand side of Eq. 4 de- 
scribes fully contained events, the second term covers lep- 
tons originating outside the detector. It may be neglected for 
electrons and tauons. 

Finally, it ought to be noted that the decay of tauons results 
in the production of both electrons and muons and hence con- 
tributes to the respective event rates. 



4 Results 

The formalism outlined in the previous section can be used 
to obtain the event rates corresponding to various of the flux 



bounds discussed in Sect. 2. This has been carried out for the 
Waxman & Bahcall and the MPR flux bounds. The results 
(for a lower energy cutoff of 100 GeV) are given in Fig. 3. 

As in Fig. 3 different bins correspond to the same solid 
angle, the form of the curves is a direct measure of attenua- 
tion effects due to crossing the Earth. This suggests that one 
might use the event rate as a function of the nadir angle in 
order to perform a tomography of the inner Earth (Jain et al., 
1999). 

To this end, we start by noting that the column density X 
may be considered as the one-dimensional Radon transform 
of the number density p: 



J p(r)S(p 



£r)d : 



r. 



(5) 



Here, p and £ are the distance from the middle of the Earth 
and the normal vector of the line of integration, respectively. 
Now, using the inversion formula for two-dimensional Radon 
transforms (Deans, 1993) one can show that X(p, £) = X(p) 
and that one may thus invert Eq. 5 by means of 

p(r) = -t—^ / X(V^T^)dq, (6) 

2?rr dr J_j R i_ r i 

where p(r) = for r > R. However, if the initial neu- 
trino spectrum (before crossing the Earth) is assumed to be 
known, one may deduce the column density directly from the 
observed event rate. Then Eq. 6 yields the density, so that a 
tomography is indeed possible. 

Unfortunately, there is an important caveat: Due to the 
statistical error a in the detection rates, the column densities 
can be known only with a rather limited precision. We take 
this into account by either adding or subtracting a throughout 
and comparing the results thus obtained with that for a = 0. 
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r/km r/km 

Fig. 4. Left: Earth density p as a function of the radius r obtained by performing a tomography on the event rates N of Fig. 3 for the 
MPR flux bound with high neutron source opacity. Shown are the results for no statistical error (solid), for iV M ± \J 7V M (dashed), and for 

iV T ± y/Tir (dot-dashed lines). A lower energy cutoff of 100 GeV is assumed. For comparison the "correct" density according to Dziewonski 
and Anderson (1981) is included (dotted line). Right: Same as left figure, but for the Waxman & Bahcall flux bound. 



The results of this analysis for a one-year run are shown in 
Fig. 4 for an MPR and for the Waxman & Bahcall flux bound. 
One can see that the maximum event rates to be expected for 
next-generation neutrino telescopes put severe constraints on 
a tomography of the inner Earth. 
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